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ABSTRACT 

Recent follow-up observations of the December 27 giant flare of SGR 1806-20 have detected a multiple- 
frequency radio afterglow from 240 MHz to 8.46 GHz, extending in time from a week to about a month 
after the flare. The angular size of the source was also measured for the first time. Here we show 
that this radio afterglow gives the first piece of clear evidence that an energetic blast wave sweeps up 
its surrounding medium and produces a synchrotron afterglow, the same mechanism as established for 
gamma-ray burst afterglows. The optical afterglow is expected to be intrinsically as bright as ma ~ 13 
at t < 0.1 days after the flare, but very heavy extinction makes the detection difficult because of the low 
galactic latitude of the source. Rapid infrared follow-up observations to giant flares are therefore crucial 
for the low-latitude SGRs, while for high-latitude SGRs (e.g. SGR 0526-66), rapid follow-ups should 
result in identification of their possible optical afterglows. Rapid multi-wavelength follow-ups will also 
provide more detailed information of the early evolution of a fireball as well as its composition. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — stars: individual (SGR 1806-20) 



1. INTRODUCTION 

Soft 7-ray repeaters (SGRs) distinguish themselves from 
classical gamma-ray bursts (GRBs) by their repetitive, 
soft bursts coming from nearby sources, likely strongly 
magnetized neutron stars dubbed "magnetars" (Thomp- 
son & Duncan 1995). Giant flares are rare events from 
SGRs, each characterized by an initial hard spike with 
more than a million Eddington luminosity and a pulsat- 
ing tail persisting for hundreds of seconds (Mazets et al. 
1979; Hurley et al. 1999; Borkowski et al. 2004). A giant 
flare was recorded on 2004 December 27 from SGR 1806- 
20 (Borkowski et al. 2004; Hurley et al. 2004; Mazets 
et al. 2004; Palmer et al. 2004; Smith et al. 2005). In 
comparison with previous ones, this flare is exceptionally 
strong with an isotropic energy of gamma-ray emission 
larger than 8 x 10 45 ergs (Boggs et al. 2005), given a source 
distance d = 15.1kpc (Corbel & Eikenberry 2004). 

Afterglow emission from giant flares was first detected 
from SGR 1900+14 at radio frequencies (Frail et al. 1999). 
The radio afterglow of the December 27 giant flare from 
SGR 1806-20 was detected about 7 days after the flare 
with a flux density at 8.46 GHz approximately 100 times 
brighter than the peak flux density of the SGR 1900+14 ra- 
dio afterglow (Gaensler et al. 2005; Cameron et al. 2005). 
Generally, the multi-frequency radio afterglow light curves 
display a power law decay behavior with time during the 
observational time from one week to one month after the 
flare, which is quite similar to the afterglows of cosmo- 
logical GRBs. But the decay rate at higher frequencies 
(e.g. 8.46 GHz) is larger than that at lower frequencies 
( e.g. 240 MHz) and the spectra clearly show a gradual 
steepening with time (Cameron et al. 2005). 

Cheng k. Wang (2003) suggested the possibility that 
the radio afterglow of the SGR 1900+14 giant flare might 



be explained with the same mechanism as established for 
GRB afterglows (Meszaros & Rees 1997), i.e. the emis- 
sion from the expanding blast wave that forms when the 
outflow from the SGR interacts with the surrounding in- 
terstellar medium (ISM). However, the data of that radio 
afterglow are too sparse for us to make a firm conclusion. 
The radio afterglow of the December 27 flare has much 
more abundant data of the light curves which was detected 
from 240 MHz to 8.46 GHz in frequency and from a week 
to a month after the flare in time. The angular size of the 
source was also measured for the first time. In this Let- 
ter, we will show that this radio afterglow gives the first 
piece of clear evidence that an energetic blast wave sweeps 
up its surrounding medium and produces a synchrotron 
afterglow. 

2. THE BLAST WAVE MODEL FOR THE RADIO 
AFTERGLOW 

It was suggested by Huang et al. (1998) that relativistic 
fireballs similar to those in classic GRBs may exist in SGR 
bursts. Thompson & Duncan (2001) proposed that the 
extraordinarily high peak luminosity L > 10 6 LEdd (where 
LEdd is the Eddington luminosity), hard spectrum and 
rapid variability of the initial spike emission of giant flares 
imply that the spike emission must originate from a rela- 
tivistically expanding fireball with the initial Lorcntz fac- 
tor of at least several in order to avoid the pair production 
problem. Such a relativistic outflow, after emitting the 
hard spike, should retain some energy and then drives a 
blast wave expanding into the surrounding medium. 

2.1. The angular size and constraint on the flare energy 

Now let's first consider the dynamics of a blast wave 
that is driven by the relativistic outflow and expands 
into the ISM. The blast wave will be gradually deceler- 
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ated by the swept-up ISM and enter the non-relativistic 
phase at t nr = (SE/ATmrripC 5 ) 1 / 3 = ASE^n^ 1 ^ 3 d, where 
E is the isotropic kinetic energy of the blast wave, c 
is the speed of light, n is the number density of ISM, 
m p is the proton mass, and we have used the usual no- 
tation a = 10 n a n in c.g.s. units except that the ob- 
server's time t is in unit of day. So, for E^/tiq < 4, 
the blast wave must have entered the non-relativistic 
phase during the observational time frame (t > 7d). 
During this non-relativistic phase, the speed (in units 
of c) and the radius of the blast wave are given by 

[3 = [12E / {UbTmmpC 5 ?)] 1 ^ = 0AE^\ 3/5 n^ 1/5 and 

R = 5[3ct/2 = 2.6 x 10 16 £^ 5 ii /5 n~ 1/5 cm respectively 
(Cheng & Wang 2003). The measured angular size at 
8.46 GHz of the radio afterglow is ~ 75 mas around seven 
to ten days after the flare and has no noticeable varia- 
tion in diameter during this period (Gaensler et al. 2005; 
Cameron et al. 2005). At a distance of d — 15.1kpc, this 
angular size implies a subluminal speed of ~ 0.35c for the 
averaged projected expansion and an apparent radius of 
R± = 9 x 10 15 cm. So, if the outflow is spherical at this time 
(i.e. R = R±), we obtain the constraint on the outflow en- 
ergy, i.e. E ~ 10 44 noergs 1 . If, instead, the outflow is 
mildly beamed in the non-rclativistic phase, then the half 
opening angle 6 satisfies sin6* = R±/R = 0.4(£46/no) -1 / 5 . 
In any case, the real kinetic energy (after beaming correc- 
tion) of the outflow is about 10 44 — 10 45 ergs for a typical 
ISM with n = lcm~ 3 and d = 15.1kpc. It is more reason- 
able to believe that the outflow becomes approximately 
spherical at the non-relativistic phase when the sideways 
expansion takes place (Livio & Waxman 2000). In the 
following we will focus on the spherical outflow case with 
E = 10 44 no ergs. The apparent radius should increase 
with time as R oc t 2 / 5 in the non-relativistic phase, so it 
will increase by a factor of 1.15 from 7 to 10 days. This 
is consistent with the fact that there is no noticeable ob- 
served variation of the angular size during the period 2 . 

2.2. The light curves of the radio afterglow 

As the blast wave sweeps up the surrounding medium, 
the shock accelerates electrons and amplifies the magnetic 
field in the swept-up matter (Blandford & Eichler 1987). 
The afterglow emission arises from synchrotron radiation 
of these shocked electrons in the magnetic field. The elec- 
tron energy distribution in the downstream of the shock 
usually manifests as a single power law or broken power 
law. The radio-band spectra show a gradual steepening 
with time, with a single power-law index of (3\ = 0.62±0.02 
around t = 7d, while at t = 23d, (3 2 = -0.91 ± 0.08 
(Cameron et al. 2005). There are also reports that a 
possible break exists in the radio-band spectrum at early 
epoches (Gaensler et al. 2005). For example, at about 9 
days after the giant flare, the measured power law spectral 
index between 4.8 and 8.6 GHz is ~ —1.0, while between 
2.4 and 4.8 GHz the index is <~ —0.66 (Gaensler et al. 
2005). The steepening of the light curve decay with fre- 
quency also favors that there is a spectral break between 



240 MHz and 8.46 GHz in the synchrotron blast wave 
model. This break is unlikely to be the cooling break be- 
cause the cooling frequency is much higher than the radio 
band (see Eq.2). As suggested by Li & Chevalier (2001) to 
explain the afterglows of GRB991208 and GRB000301C, 
we consider that this break is caused by an intrinsic break 
in the energy distribution of the shock-accelerated elec- 
trons, which takes the form of 
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where 7 m i n and 75 are the minimum and break Lorcntz 
factors respectively. The observed synchrotron spectrum 
of the Crab Nebula supports the intrinsic break in the in- 
jected particle spectrum (Amato et al. 2000). Similarly, 
we assume that the ratio of the break energy and the mini- 
mum energy, Rb = ib/lmim remains essentially a constant 
in time so that the relative shape of the distribution is 
time invariant (Li & Chevalier 2001). From the values of 
[3\ and /?2, we infer p\ ~ 2.2 and p 2 — 3.0. 

Assuming that shocked electrons and the magnetic field 
acquire constant fractions (e e and es) of the total shock 
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energy density, we get 7 min 

(16ne B f3 2 nm p c 2 ) 1 / 2 = 1.4 
Thus we can obtain the break frequency corresponding to 
7m i n , the break frequency corresponding to the electron 
break caused by the synchrotron cooling, and the peak 
flux of the spectrum, which are given by 
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F Vm = 95< 5 ^ /5 ^ /10 e^ /2 _ 1 (d/15.1kpc)- 2 Jy. 



For a broken power-law electron distribution described 
by Eq.(l), the light curves at the frequencies below and 
above the break frequency Vf, (which is the synchrotron 
frequency corresponding to 7&) are different. According to 
the afterglow theory, the light curves in the non-relativistic 
phase should have the decay slopes of (21 — 15pi)/10 and 
(21 — 15p2)/10, respectively. For p\ = 2.2 and p 2 = 3, the 
predicted slopes are in agreement with the observed ones 
(i.e. ai = -1.43±0.17for240MHzanda 2 = -2.63±0.13 
for 8.46 GHz) of the radio afterglow data ( Cameron et al. 
2005), if the break frequency locates between 240 MHz 
and 8.46 GHz. The 4.86 GHz light curve shows a break 
around t — 9d, which implies that, in the above picture, 
Vb just crosses 4.86 GHz at this time. This then suggests 
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The shifting of Vb to 



lower frequencies would make 1.43 GHz light curve decay 
faster at late times. 

The flux density at 1.43 GHz is 
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lr The discrepancy between the isotropic gamma-ray energy and the kinetic energy of the afterglow may be naturally expected within the 
magnetar framework as the prompt hard spike emission taps the magnetic energy, while the afterglow is only related to the kinetic energy 
(Zhang & Kobayashi 2005). On the other hand, this discrepancy could be alleviated by a denser environment. 

2 After the submission of the paper, more observation data of the angular size at time from 10 to 20 days are reported in Gaensler et al. 
(2005). The angular size shows a slow increase with time, being consistent with the _R oc t 2 / 5 relation. 
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where the coefficient is calculated for p\ = 2.2. At t = 
9.87d, the observed flux density at 1.43 GHz is 96 ± 2mJy 
(Cameron et al. 2005). So, for typical parameter values of 
E = 10 44 ergs, n = 1cm -3 , d = 15.1kpc and e e = 0.3, we 
infer e# = 0.04. This equipartition factor is close to that 
inferred for GRB afterglows (e.g. Wijers & Galama 1999; 
Panaitescu & Kumar 2001) and corresponds to a magnetic 
field of ~ 10 mG in the shocked medium. Nevertheless, we 
should note that the observations are all in the power-law 
domain and the spectral transitions (e.g. v m and v c ) have 
not been observed, which makes the physical parameters 
of this flare not well constrained. 

We have also carried out detailed numerical calculations 
of the radio afterglow light curves at five frequencies, i.e. 
8.46 GHz, 4.86 GHz, 1.43 GHz, 608 MHz and 240 MHz. 
The numerical model contains more careful treatments to 
the blast wave dynamics during the trans-relativistic stage 
(Huang et al. 2000) and the broken power-law electron 
distribution 3 . It is found that the numerical model can 
give a satisfactory fit to the light curve data over five radio 
frequencies, as shown in Fig.l. The theoretic light curve of 
610MHz are somewhat above the observation data. This 
discrepancy is also evident in the power-law fit of the ra- 
dio spectra in Fig.l of Cameron et al. (2005). The pa- 
rameter values used in the fit are consistent with those 
constrained by the above analytic study 4 . An interesting 
finding is that, unlike GRB radio afterglows that usually 
peak around a week after the burst, the predicted SGR 
radio afterglows are significantly stronger at earlier times. 
This calls for an even more rapid follow-up in the radio 
band. 

3. OPTICAL, INFRARED AND X-RAY AFTERGLOWS OF 
GIANT FLARES AND INTERMEDIATE BURSTS OF SGRS 

SGR afterglows have been detected only in the radio 
band so far. Within the above picture, it is natural to 
expect optical and infrared afterglows as well produced by 
the same relativistic fireball accompanying the SGR giant 
flares. At the early time t < i nr , the blast wave is rela- 
tivistic and its emission spectrum is characterized by 
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The optical flux density at -ft-band (4.55 x 10 Hz) is 
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spherical phase, the above estimate is conservative. Be- 
cause the characteristic synchrotron frequency v m at early 
times is generally below the optical frequency for an initial 
Lorentz factor of T <~ 10, the optical afterglow reaches 
its peak at the deceleration time of the relativistic out- 



flow, which is tde 
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afterglow has a flux density of 
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where p ~ 2.4 is the power-law index of the electron distri- 
bution inferred from the global spectrum F v oc j,- - 68 ± -°2 
(Cameron et al. 2005). For typical parameter values and 
the inferred value cb = 0.04 from the radio afterglows, 
the magnitude of the optical afterglow at t = O.ld is 
about tur ~ 13. Considering that, for an initial beamed 
outflow, the isotropic equivalent energy in the relativis- 
tic phase may be larger than that in the non-relativistic 

3 The code also takes account of the re-derived electron distribution by 

4 For a smaller source distance d (Cameron et al. 2005) , E/n should be 
angular size. As a result, a larger value for e_g is then inferred. However, 
wave model remains unchanged. 



Such a bright optical afterglow would be easily detected 
by rapid response detectors, such as Ultraviolet/Optical 
Telescope (UVOT) on Swift and Robotic Optical Tran- 
sient Search Experiment (ROTSE), if the optical extinc- 
tion to the source is not very large. However, SGR 1806-20 
has a low Galactic latitude and has optical extinction of 
Ay <~ 30 mag (Eikcnberry et al. 2001). This makes it im- 
possible to observe the optical afterglow from this SGR. 
The difficulty also applies for the other two SGRs at low 
Galactic latitudes, SGR 1900+14 and SGR 1627-41, whose 
extinction is also large (Vrba et al. 2000; Wachter et al. 
2004). The i?-band extinction of SGR 1900+14 is moder- 
ate , i.e. A ft ~ 7 mag (Akerlof et al. 2000), leaving some 
hope for optical afterglow detection at its peak. The best 
target for optical follow-up observations among the four 
known SGRs is SGR 0526-66, which lies in the Large Mag- 
ellanic Cloud and has Ay — 1. At a distance of d ~ 50 kpc, 
the optical afterglow of a similar giant flare from this SGR 
is expected to be bright enough for follow-ups. Even the 
less rare intermediate bursts with the isotropic energies 
ranging from 10 41 to 10 43 ergs (e.g. Golenetskii et al. 
1984; Guidorzi et al. 2004) are expected to have optical 
afterglows brighter than rriR ~ 23 at t < O.ld for SGR 
0526-66. This is well above the Swift UVOT sensitivity. 
Detecting such an optical afterglow would further test the 
blast wave model and provide more detailed information 
of the early evolution of a fireball. 

Because of the high optical extinction of low-latitude 
SGRs, infrared observations to giant flares and interme- 
diate bursts appear very useful. The K band flux den- 
sity would be, according to Eq. (6), ^2.2/im ~ 3mJy or 
rriK — 13 at t = Id for typical parameter values to model 
the radio afterglow of the December 27 giant flare. Such a 
bright infrared afterglow could be well detected by ground- 
based infrared telescopes. An example of the (unabsorbed) 
light curves of the infrared and optical afterglows calcu- 
lated with the numerical model for the radio afterglows 
is also shown in Fig. 2. The flux of the X-ray afterglow 
can be estimated in the same way. For the above typical 
parameter values, we find the flux at t — O.ld is about 
F\ ~ 2.6 x 10 -10 ergcm _2 s _1 . This flux is, however, much 
dimmer than the observed X-ray flux of the SGR 1900+14 
giant flare (Woods et al. 2001), which is thought to be 
originated from the neutron star surface immediately af- 
ter the burst. 

Huang & Cheng (2003) for afterglows in the deep Newtonian phase. 

correspondingly smaller in order to be consistent with the measured 
the global agreement of the radio afterglow behavior with the blast 
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4. DISCUSSIONS 

We have shown that the radio afterglow behavior of the 
December 27 giant flare of SGR 1806-20, including the an- 
gular size and multi-frequency light curves, provides a first 
piece of clear evidence that the fireball/blast wave model 
indeed works in SGRs as well. The polarization emis- 
sion detected in the radio afterglow (Gaensler et al. 2005; 
Cameron et al. 2005) at a comparable level (1% — 3%) to 
GRB afterglows (Covino et al. 1999; Wijers et al. 1999) 
is also consistent with the blast wave model. The initial 
Lorentz factor of the outflow that drives the blast wave 
is not well constrained, because only the late-time radio 
afterglow is available up to now. The composite of the 



relativistic outflow is also unknown and could be normal 
baryon matter as well as relativistic electron-positron dom- 
inated wind since both could drive relativistic blast waves 
when they are interacting with the surrounding medium. 
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days since the SGR 1806-20 giant flare 



Fig. 1. — Numerical fits of the radio afterglows of the December 27 giant flare from SGR 1806-20 with the blast wave model. The data 
are taken from Gaensler et al. (2005) (open symbols) and Cameron et al. (2005) (solid symbols). The parameter values used in the fits are 
E = 10 44 erg, n = 1cm -3 , T = 10, e e = 0.32, e B = 0.1, pi = 2.2, p 2 = 3.2, R b = 150 and d = 15.1kpc. The inset shows these light curves 
over a longer time range. 
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Fig. 2. — The predicted optical (i2-band) and infrared (J^-band) afterglows of a giant flare similar to the December 27 event. Parameter 
values used are the same as those in the fits of Fig. 1. 



